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ABSTRACT 

Because accretion and merger shocks in clusters of galaxies may accelerate particles to high energies, 
clusters are candidate sites for the origin of ultra-high-energy (UHE) cosmic-rays. A prediction was 
presented for gamma-ray emission from a cluster of galaxies at a detectable level with the current 
generation of imaging atmospheric Cherenkov telescopes. The gamma-ray emission was produced via 
inverse Compton upscattering of cosmic microwave background (CMB) photons by electron-positron 
pairs generated by collisions of UHE cosmic rays in the cluster. We observed two clusters of galaxies, 
Abell 3667 and Abell 4038, searching for very-high-energy gamma-ray emission with the CANGAROO- 
HI atmospheric Cherenkov telescope system in 2006. The analysis showed no significant excess around 
these clusters, yielding upper limits on the gamma-ray emission. From a comparison of the upper limit 
for the north-west radio relic region of Abell 3667 with a model prediction, we derive a lower limit 
for the magnetic field of the region of ~0.1/iG. This shows the potential of gamma-ray observations 
in studies of the cluster environment. We also discuss the flux upper limit from cluster center regions 
using a model of gamma-ray emission from neutral pious produced in hadronic collisions of cosmic- 
ray protons with the intra-cluster medium (ICM). The derived upper limit of the cosmic-ray energy 
density within this framework is an order of magnitude higher than that of our Galaxy. 
Subject headings: galaxies: clusters: individual(Abell 3667, Abell 4038) — gamma rays: observations 
— galaxies: magnetic fields 
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Clusters of galaxies are the largest systems in the Uni- 
verse that are gravitationally bound, and they are po- 
tential sources of ultra-high-energy (UHE) cosmic rays, 
since their sizes and moderate magnetic fields allow a 
high maximum energy (^10^'' eV) in acceleration (Os- 
trowski 2002). Although cluster accretion and merger 
shocks could produce such high-energy particles, accre- 
tion shocks may be more effective than merger shocks 
in particle acceleration, due to their high Mach num- 
bers (Miniati et al. 2000; Ryu et al. 2003). Cosmic- 
ray electrons accelerated directly by these shocks may 
produce gamma-ray emission via inverse Compton (IC) 
scattering of the cosmic microwave background (CMB) 
(Totani & Kitayama 2000; Miniati 2003; Gabici & Blasi 
2004). On the other hand, accelerated cosmic ray pro- 
tons can interact hadronically with the intra-cluster 
medium (ICM), and gamma-rays may be produced via 
7r°-decay (Volk et al. 1996; Berezinsky et al. 1997; Pfrom- 
mer & Enfilin 2004) as well as IC emission by sec- 
ondary electron/positron pairs from 7r^-decay (Blasi & 
Colafrancesco 1999). 

Observations of clusters of galaxies at various wave- 
lengths (e.g., radio, EUV, X-ray) suggest the existence 
of non-thermal particles in these gigantic objects (Fusco- 
Femiano et al. 2001; Nevalainen et al. 2004; Bowyer 
et al. 2004; Giovannini & Feretti 2004). However, at 
gamma-ray energies, no observational evidence has been 
reported from clusters of galaxies (Reimer et al. 2003), 
though there is suggestive evidence (Kawasaki & Totani 
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2002; Scharf & Mukhcrjcc 2002). Observations in the 
TcV energy band with Imaging Atmospheric Cherenkov 
Telescope (hereafter lACT) experiments have yielded the 
only upper limits to date (Hattori et al. 2003; Fegan et 
al. 2005; Perkins et al. 2006; Domainko et al. 2007). 

Recently, Inoue et al. (2005) considered protons, ac- 
celerated up to 10^*^ 10^^ cV by accretion shocks around 
a massive cluster, interacting with the CMB photons, 
with secondary electron-positron pairs produced in the 
p-7 process, boosting those photons into the TeV en- 
ergy range by IC scattering. Although their predic- 
tion depends on many physical parameters, the pre- 
dicted gamma-ray flux could be at a detectable level 
for current lACT experiments for massive, nearby clus- 
ters. Thus, observations of clusters of galaxies with 
lACTs probe high-energy processes and the environment 
in these large-scale systems, and if gamma-ray signals arc 
detected, they may also provide clues to help solve the 
mystery of UHE cosmic-ray production. If there is no 
detected signal, we can place limits on the physical pa- 
rameters of clusters, such as the strength of the magnetic 
field. 

In this paper we report on a search for TeV gamma- 
ray emission from two clusters of galaxies, Abell 3667 
and Abell 4038, with CANGAROO-III, an array of imag- 
ing atmospheric Cherenkov telescopes. We selected these 
targets according to their high masses and relative close- 
ness from the southern Abell catalog (Abell et al. 1989). 

Abell 3667, also known as SC 2009-57, is classified as 
type L in the Rood-Sastry system (Rood & Sastry 1971) 
due to the linear arrangement of the galaxies, including 
two of the brightest D galaxies. The cluster has a red- 
shift oi z — 0.055 (Sodre et al. 1992), and is centered at 
[a(2000) = 20''12™27".4, (5(2000) = -56°49'36"], which 
is the location of the brightest D galaxy (Knopp et al. 
1996). The cluster is one of the brightest X-ray sources in 
the southern sky (Edge et al. 1992), and is also known to 
show significant diffuse radio emission around its center 
(Rottgering et al. 1997). 

Abell 4038, also known as Klemola 44, is a rich south- 
ern cluster with z = 0.028, and is classified as type cD 
in the Rood-Sastry system, where the cD galaxy is cen- 
tered at [a(2000) = 23''47™45M, (5(2000) = -28°08'26"] 
(Slee et al. 2001). An X-ray image shows an extcndeci 
morphology, and there is a radio relic near the cD galaxy 
(Slee and Roy 1998), though it is smaller than the point 
spread function of typical lACTs (^0.1°). 

The search described in the following sections focused 
on the detection of point sources within the cluster fields 
as well as looking for gamma-ray signals from several 
regions by assuming gamma-ray emission models: The 
giant radio relics around Abell 3667 may indicate the 
sites of shocks, where particle acceleration occurs effec- 
tively, and thus ganuua-ray emission could be expected 
from the relics by the scenario of the UHE proton ori- 
gin (Inoue et al. 2005). Apart from the shock regions, 
the density of the ICM is highest at the cluster centers, 
so the gamma-ray flux via Tr'^-decay would be strongest 
there. 

In §2, we introduce the CANGAROO-III telescope 
systems and observations of the clusters. The data- 
analysis procedures are explained in §3, and the main 
results together with a definition of the gamma-ray 



search regions are described in §4. Finally, a discus- 
sion of the gamma-ray emission from clusters of galax- 
ies based on the CANGAROO-III results is presented. 
Throughout this paper, we assume a Hubble constant of 
Ho=70kms-iMpc-^ 

2. INSTRUMENT AND OBSERVATIONS 

Two clusters of galaxies, Abell 3667 and Abell 
4038, were observed in the TeV energy band using 
the imaging atmospheric Cherenkov technique with 
the CANGAROO-III telescope system (Enomoto et al. 
2006a), located near Woomcra, South Australia (136.786 
degree E, 31.099 degree S, 160m a.s.L). The system con- 
sists of four telescopes, which are located at the corners of 
a diamond with 100 m sides (Enomoto et al. 2002). The 
first telescope, which we call Tl, was not used in these 
observations, since its current performance is inferior to 
that of the other telescopes. The specifications of the 
second, third, and fourth telescopes, hereafter called T2, 
T3, and T4, are almost the same: they have segmented 
paraboloid reflector 10 m in diameter and 8 m in focal 
length. Each reflector is composed of 114 spherical mir- 
ror facets, and the total mirror area is 54 m^ (Kawachi 
et al. 2001). At each focus there is an imaging camera, 
which is an array of 427 photomultipliers (PMTs). Each 
PMT covers a sky field of 0.17° in diameter, and the total 
field-of-view is about 4°, suitable for applying analysis 
using an imaging technique (Kabuki et al. 2003). The 
data acquisition system is triggered when at least two 
telescopes have signals coinciding for more than lOnsec 
within a 650nsec time window, thus eliminating muon 
events that mostly trigger a single telescope (Nishijima 
et al. 2005). Then, the amplitude and the arrival times of 
signals from PMTs are digitized by ADC/TDC modules, 
and recorded for off-line analysis (Kubo et al. 2001). 

The observations were carried out for ON-source and 
OFF- source tracking nms. For OFF-sourcc nms, the 
target position was shifted in right ascension so that the 
telescopes tracked the same trajectory across the sky as 
ON-source runs. Also, we adopted wobble mode observa- 
tions for both ON-source and OFF-source runs, in which 
the pointing direction was shifted in declination by ±0.5° 
from the target direction every 20 minutes. One of the 
advantages of the wobble mode is to average the response 
of PMTs, since the target position rotates on the FOV. 
All observations were made on moonless nights from July 
to September, 2006. Details of the data sets are summa- 
rized in Table 1. 

In addition to the observations, we also observed dark 
regions (with no bright stars or gamma-ray sources in 
the field-of-view) without imposing the telescope coinci- 
dence described above, and we extracted local muon-ring 
events from these data to monitor the total performance 
of the telescopes (Enomoto ct al. 2006a). This calibra- 
tion (denoted muon run) was done every month, and 
their statistics are also shown in Table 1. 

3. DATA ANALYSIS 

We followed the analysis procedure explain(;d in detail 
in Enomoto et al. (2006a), which we briefly describe here. 

First, raw data were calibrated, using daily calibra- 
tion runs using LEDs (Kabuki ct al. 2003). We then 
selected the shower events from the calibrated datasets. 
For each triggered event, camera pixels that recorded less 
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TABLE 1 

Summary of the data set used in this analysis. 



Cluster 


Term 


ON'' 


OFF'' 








[hour] 


[hour] 


[deg] 


Abell 3667 


Jul- Aug (2006) 


29.7 


23.7 


28.6 


Abell 4038 


Aug-Sep (2006) 


23.6 


17.7 


13.1 


muon run 


Jul-Sep (2006) 


22.4 




10.8 



'''Dcad-tiiric corrected observation time used in our analysis. 
''Mean zenith angle of ON-source runs. 

than 5 p.e. were discarded so as to remove any night-sky- 
background photons, and shower events were extracted 
from the remaining pixels, imposing the conditions that 
there were at least 5 adjacent hit pixels, and that the 
pixels were triggered within ±30 nsec from the mean ar- 
rival time of all hit pixels. This procedure cleaned the 
shower images and separated random noise, such as mul- 
tiple night-sky-background photons. The typical shower 
event rate was ^7 Hz (average over 5 minutes), and wc 
excluded the data from the analysis when the shower 
event rate was below 5 Hz so as to remove any data af- 
fected by clouds etc. The effective total observation times 
for the selected datasets, taking account of the dead time 
in data acquisition, are summarized in Table 1 . After this 
image-cleaning procedure, wc discarded events with any 
hits in the outermost layer of the imaging cameras since 
such shower images may be distorted (Enomoto et al. 
2006b). 

Next, image moments (width and length) of showers 
were calculated as defined by Hillas (1985), and the ar- 
rival direction of the shower was reconstructed event by 
event, by minimizing the sum of the squared widths of 
the images weighted by their total photo-electron num- 
bers seen from the assumed direction, as described in 
Kabuki et al. (2007). 

Finally, gamma-ray/hadron separation was carried out 
by applying the Fisher discriminant method (Fisher 
1936; Enomoto et al. 2006a). In this method, the Fisher 
discriminant (hereafter FD) is defined as a linear combi- 
nation of the image moments, 

6 

FD=J2ai-Pi, (1) 

where P = (W2, W3, W4, Z/2, ^3, -^^4) is a set of energy- 
corrected width and length of shower images of three 
telescopes. The coefficients, ai [i — 1^6), were deter- 
mined so that the; difference of the FD distribution of 
gamma-ray events and that of hadron events would be 
maximized. We used Monte-Carlo simulation data as 
gamma-ray events, and OFF-source run data as back- 
ground hadron events for deriving the coefficients. We 
then extracted gamma-ray events from ON-source run 
data by fitting the ON-source FD distribution with that 
of the background (OFF-source) distribution plus the 
gamma-ray distribution. Our Monte-Carlo simulation 
code is based on GEANT3, the details are described in 
Enomoto et al. (2002), where some parameters such as 
the geometry of the telescopes are replaced with those 
of the current CANGAROO-HI system. The degrada- 
tion of the overall light collection efficiency (including 
refiectivities of the refiectors, quantum efficiencies of pho- 
tomultipliers, etc.) and the spot size of each telescope 



were estimated from a muon ring analysis (Enomoto et 
al. 2006a) using muon run data, and they were included 
in our simulation. In the gamma-ray simulation, we as- 
sumed a power-law spectrum index of 7 = —2.1, which 
is often assumed for clusters of galaxies (e.g., Volk & 
Atoyan (2000)). 

4. RESULTS 

4.1. Two dimensional morphology and 6^ distribution 

First, we calculated the two-dimensional (2D) signifi- 
cance map around the cluster centers. We divided these 
(ON source) regions into 0.2° xO. 2° square bins, and 
calculated the gamma-ray like excesses and their errors 
with the FD fitting method, described in the previous 
section. Each background (OFF source run) bin was 
taken so that its position on the field of view would cor- 
respond to that of the ON region's bin, but the area 
was extended to 3x3 neighboring bins, to improve the 
statistical accuracy. Fig. 1 shows the resulting 2D sig- 
nificance maps of gamma-ray like excesses. Since the 
gamma-ray acceptance falls off toward the outer part of 
the field-of-view, we limit the map to within 1 degree 
from the cluster centers. The significance distributions 
from all bins in 2D maps were well approximated by stan- 
dard normal distributions for both regions. The best fit 
Gaussians have mean values of 0.02±0.10 (Abell 3667) 
and 0.17±0.11 (Abell 4038), with standard deviations of 
1.08±0.08 (Abell 3667) and 1.02±0.08 (Abefi 4038), and 
there are no significant gamma-ray signals. 

Next, we show the 6*^ distribution, where 9 is the space 
angle between the target position and the reconstructed 
arrival direction, from the cluster centers in Fig. 2. The 
background bin for calculating each 6^ bin was taken 
from the OFF source region to correspond to the po- 
sition in the field of view. Although there were devia- 
tions in the 9"^ distributions, they were not significant 
(<3(t), considering our point-spread function {9^ < 0.06 
degree^). 

In summary, there were no detectable point sources in 
the cluster fields. 

4.2. Gamma-ray emission profiles and upper limits 

We also adopted several gamma-ray emission profiles 
in the cluster fields, and searched for diffuse gamma-ray 
emission. First, we defined two circular regions (here- 
after NW/SE Relic regions) that cover the prominent 
radio relics around Abell 3667, since they may represent 
a shock morphology. The center coordinates (R.A. & 
Dec. in J2000) and their radii were defined as follows: 
(20''10'"24^ -56°27'00") and 0.30° for the NW Relic re- 
gion, (20''14™36% -57°03'00") and 0.24° for the SE Relic 
region. 

We expect that gamma-ray emission via 7r°-decay is 
concentrated at the cluster center regions. It is well 

known that many clusters have a diffuse X-ray morphol- 
ogy at their centers, which may trace thermal compo- 
nents bounded by the gravitational potential of clusters. 
We thus assume that the gamma-ray emission profile 
traces the X-ray morphology of clusters. We adopted 
the ROSAT PSPC data for the X-ray morphology The 
peak positions of the X-ray brightness are almost coin- 
cident with the cD galaxies of the clusters, which were 
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Fig. 1. — Two-dimensional significance maps around the clusters. The map centers correspond to the position of the cD galaxy of 
each cluster, and at each point a smoothing algorithm within 0.6° xO. 6° was applied. Contours of other wavelengths (X-ray & radio) are 
over-plotted. The left panel is for Abell 3667, where the white contours show ROSAT hard-band data (Voges et al. 1999) and orange 
contours show SUMSS 843 MHz radio data (Mauch et al. 2003). The right panel is for Abell 4038, where white contours show ROSAT 
hard-band data and orange contours show VLA 1.4 GHz radio data (Condon et al. 1998). Our point-spread function is also shown at the 
bottom right-hand corner. 

which can be regarded as being a point source, consid- 
ering the positional resolution of CANGAROO-III. So 
we first searched for gamma-ray signals from the Abell 
3667 center region based on ROSAT data, and we later 
discuss the case of a point source, especially concern- 
ing the cosmic ray energy density. With the Tr'^-decay 
model, Volk & Atoyan (2000) assumed that the high- 
energy protons were accumulated in a cluster through su- 
pernova explosions, and that the predicted proton spec- 
trum forms the power-law index F — —2.1 with an energy 
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Fig. 2. — 9^ distributions from cluster centers (where the def- 
inition of the center is the same as in Fig. 1). The upper panel 
is for Abell 3667 and the lower panel is for Abell 4038. The solid 
histograms show Monte-Carlo simulation results, assuming a Crab- 
nebula-level point-like gamma-ray source at each cluster center, 
and the dotted lines indicate the zero point. Note that our PSF 
corresponds to 9^ < 0.06 degree^ . 

the tracking points of our observations. We then defined 
two regions (hereafter Cluster Core regions), such that 
their centers were at the position of the cD galaxies and 
the radii were equal to the point where the S/N of the 
ROSAT data fell below - 3, which was 0.40° for Abell 
3667 and 0.26° for AbeU 4038, as described in Table 2 of 
Mohr et al. (1999). For Abell 3667, more recent obser- 
vations with higher resolution have been reported (e.g., 
XMM-Newton (Briel et al. 2004); Chandra (Vikhhnin et 
al. 2001)). Chandra has a limited field-of-view for our 
purpose, however XMM showed that the signal region 
above the background noise level was in the central 11' 



tion, gamma-rays were generated uniformly within the 
defined area with a power-law spectrum having an index 
of 7 = -2.1. 

The gamma-ray events from each region were calcu- 
lated by the FD fitting method, as before. The FD dis- 
tribution of each region fitted with that of OFF-source 
region and that of gamma-ray events from a Monte-Carlo 
simulation is shown in Fig. 3. A gamma-ray signal would 
appear around FD=0 (see, e.g., Fig.l in Enomoto et al. 
(2006b)); however, the calculated significances from the 
4 regions did not exceed 3a, and so there is no evidence 
of extended emission. 

Therefore, we calculated the 2a upper limits on the 
integral gamma-ray fluxes from these regions. The ob- 
tained fiux upper limits, their threshold energy, and the 
definition of each region are summarized in Table 2. 

5. DISCUSSION 

Inoue et al. (2005) predicted gamma-ray emission at 
accretion shocks around a massive cluster. We searched 
for gamma-ray emission from radio relics of Abell 3667, 
assuming that they might trace the accretion shock 
(EnBlin et al. 1998), although it has also been suggested 
that the relics are the results of a major merger (Roet- 
tiger et al. 1999) in which case the particle accelera- 
tion would not be as strong as assumed in the accretion 
model. We found no evidence of gamma-ray emission 
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TABLE 2 

Summary of results; assumed gamma-ray emission region, threshold energy, and 2cr upper limits. 



Chistcr 



Confer 



Radius Tlirosliold 2(t ti 







li. A. (.12000) 


DEC.: (.12000) 






[cm "bc'c 




Abell 6obi 


iV W Helic 


2U iU 24 


—00 27 UU 


0.30 


0.95 


3.19X 10~ 


12 












1.45 


1.64x10" 


12 












2.05 


1.05x10- 


12 




SE Relic 


20'! 14^308 


-57°03'00" 


0.24 


0.85 


5.69x10- 


12 
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1.86xl0~ 
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1.55x10- 
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Fig. 3. — FD distributions from tfic four regions. The crosses 
are FD distributions of the ON region, the dashed histograms are 
those of the OFF region, the dotted crosses arc their subtraction 
(ON— OFF), and the solid histogram is that of gamma rays by 
Montc-Carlo simulation, which were derived from the fitting pro- 
cedure : a) the Abell 3667 NW Relic region, b) the Abell 3667 SE 
Relic region, c) the Abell 3667 Cluster Core region, d) the Abell 
4038 Cluster Core region. The ON— OFF and gamma ray distribu- 
tions are magnified bv 2.5 for clarity. The resulted gamma-ray flux 
were at -13%, 13%, 13% and -0.9% of the Crab-level flux for (a)- 
(d), where the negative excesses represent statistical fluctuations, 
and none of the significances exceed 3<t. 

from either region. Since estimations of the magnetic 
field of Abell 3667 have been made for the area of the 
cluster center and the north-west relic so far (Johnston- 
HoUitt 2003), we compared the derived upper limits from 
NW Relic region with the model prediction, as shown in 
Fig. 4. The model assumes a proton luminosity of one 
tenth of the kinetic energy flux through strong accretion 
shocks, which depends on the cluster mass in the form of 
ocM^/^ (see Eq.(2) in Inoue et al. (2005)), and we scaled 
the predicted gamma-ray flux according to the mass (M) 
and distance (d) of Abell 3667 from the (Coma- like clus- 
ter) parameters used in their model (M=2xlO^^M0, 
d=100 Mpc). The mass of Abell 3667 has been estimated 
using the Virial relation to be 3.7xlO^^M0 (Sodre et al. 
1992) or 1.7xlOi^M0(Girardi et al. 1998), so here we 
adopt a cluster mass of their mean value, 2.7x10"'^^Mq. 
The scaled fluxes are shown in Fig. 4 with lines for mag- 
netic fields of 0.1/xG, 0.3/uG, and 1.0/iG. 



^ 1 — I — ^ 

B=0.1[nG] 
B=0.3[nG] 
B=1.0[nG] - 
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Photon energy [TeV] 

Fig. 4. — Derived gamma-ray flux upper limits from the NW 
Relic region of Abell 3667 (filled squares) with the predicted 
gamma-ray fluxes by Inoue et al. (2005). The model was scaled 
with the mass and the distance of Abell 3667 to that of a Coma- 
like cluster, and is shown in the case of magnetic fields of 0.1//G, 
0.3mG, and LO/iG. 

Fig. 4 indicates that we can set a lower limit for the 
magnetic field in the cluster to be ~ O.l/uG, within the 
framework of the model by Inoue et al. (2005). This value 
is not a strong constraint on the magnetic field when it is 
compared with other estimates, e.g., a few /iG from Fara- 
day rotation measurements (see Johnston-HoUitt (2003) 
for other results); however, the result provides an in- 
dependent method of a magnetic field estimation, using 
TeV gamma-ray observations. The fiux upper limits from 
the SE Relic region were higher than that of the NW 
Relic, and the lower limit of the magnetic field strength 
was estimated to be ~ 0.1/xG, depending on the assumed 
cluster mass. Note that the above flux upper limits also 
provide a constraint on the gamma-ray emission via pri- 
mary electron IC omission, which is believed to appear 
at the shocks (Miniati 2003). 

We also searched for gamma-ray emission from the 
Cluster Core regions, deriving flux upper limits. The 
gamma-ray flux via 7r°-decay, produced in hadronic col- 
lisions of high-energy protons with the ICM, is thought 
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to be brightest at the cluster centers, and its flux level 
is usually discussed from the perspective of an effec- 
tive confinement of cosmic-rays inside a cluster during 
the Hubble time. Here, we discuss the total cosmic- 
ray energy stored inside the cluster centers, using the 
CANGAROO-HI result. We plotted the flux upper lim- 
its from the Abcll 4038 Cluster Core region together with 
the EGRET upper hmit (Reimer et al. 2003) in Fig. 5. 
We adopted the assumptions of Volk & Atoyan (2000), as 
introduced in a previous section, and the gamma-ray ab- 
sorption effect by IR photons (P0.4 model in Aharonian 
et al. (2006)) was also incorporated. The gamma-ray 
spectra were represented by lines in Fig. 5, which were 
scaled to be consistent with the EGRET upper limits. 

As shown in Fig. 5, the EGRET & CANGAROO-HI 
results for Abell 4038 gave almost the same constraint on 
the gamma-ray emission for the case of F = —2.1, and the 
total cosmic-ray energy to explain our flux upper limits is 
1.2x10^'^ erg, using an ICM density of 10~^ cm"'^, which 
is a typical value for cluster centers (Blasi et al. 2007). 
We then derived the upper limit of the cosmic-ray energy 
density within the Cluster Core region of Abell 4038, as 
~40eV cm~'^, assuming a spherical symmetry with the 
radius that we defined for this region. The same esti- 
mation was applied for the Cluster Core region of Abell 
3667 as well, and the upper limit was calculated to be 
'^20eV cm~'^. Two factors should be considered for this 
value: as described in §1, the morphology of Abell 3667 
is elongated towards two radio relics, rather than a sim- 
ple sphere, so the assumption of a spherical symmetry 
might need to be reexamined where the cluster type of 
Abell 4038 is cD. Also, if we adopt the XMM results, the 
search region is effectively a point source, as described 
beforehand. In this lower flux upper limit is ob- 

tained, and the total volume of the search region also 
decreased, with the cosmic-ray energy density increasing 
to '^40eV cm~'^, which was the same level for Abcll 4038. 
In any case, the derived values are 1 order of magnitude 
higher than that of our Galaxy, ~1 eVcm~^, opening the 
door to discussions of the non-thermal component in the 
clusters. All the calculated upper limits are summarized 
in Table 3. 

Further observations of clusters of galaxies in the TeV- 
band by next-generation lACTs currently in the planning 
stage, such as CTA^^ or AGIS^'' , and in the GeV-band 
by GLAST^° will open a new window in the research of 
high-energy phenomena of clusters of galaxies with their 
improved sensitivities. 

6. CONCLUSION 

We observed two clusters of galaxies, Abell 3667 
and Abell 4038, searching for very-high-energy gamma- 
ray emission with the CANGAROO-HI atmospheric 
Cherenkov telescope system in 2006. No significant ex- 



cess was detected from the clusters, and flux upper limits 
on the gamma-ray emission were obtained. By compar- 
ing the upper limit for the north-west radio relic region 
of Abell 3667 with a model prediction, we can derive a 
lower limit for the magnetic field of the region of ~0.1/iG. 
We also discussed the fiux upper limit from the cluster 




10 10^ 10^ 
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Fig. 5. — Derived gamma-ray flux upper limits from the Clus- 
ter Core region of Abell 4038 with the gamma-ray emission spec- 
trum via TT^-decay process, normalized to the EGRET upper limits 
(Reimer et al. 2003). The EGRET upper limits are indicated by 
arrows at 100 MeV, and the lines are the gamma-ray spectrum 
from the proton power-law indices of F = —2.1 and —2.3 with an 
energy cutoff of 200 TeV. The gamma-ray absorption effect by IR 
photons are represented by dot-dash lines, where the PC. 4 model 
in Aharonian et al. (2006) is adopted. 



center regions using a gamma-ray emission model via the 
decay of 7r° produced in hadronic collisions of cosmic-ray 
protons with the ICM. The upper limit of the cosmic- 
ray energy density stored within cluster centers was es- 
timated to be ~40eV cm~^ by imposing some assump- 
tions, such as the ICM density, and the values are 1 order 
higher than that of our Galaxy. These estimations show 
the potential of gamma-ray observations in studies of the 
cluster environment. 
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